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ABSTRACT 

We consider the evolution of a supermassive black hole binary (SMBHB) 
surrounded by a retrograde accretion disk. Assuming the disk is exactly in the 
binary plane and transfers energy and angular momentum to the binary via 
direct gas accretion, we calculate the time evolution of the binary’s semi-major 
axis a and eccentricity e. Because the gas is predominantly transferred when the 
binary is at apocenter, we find the eccentricity grows rapidly while maintaining 
constant a(l + e). After accreting only a fraction of the secondary’s mass, the 
eccentricity grows to nearly unity; from then on, gravitational wave emission 
dominates the evolution, preserving constant a( 1 — e). The high-eccentricity 
waveforms redistribute the peak gravitational wave power from the nHz to /iHz 
bands, substantially affecting the signal that might be detected with pulsar timing 
arrays. We also estimate the torque coupling binaries of arbitrary eccentricity 
with obliquely aligned circumbinary disks. If the outer edge of the disk is not an 
extremely large multiple of the binary separation, retrograde accretion can drive 
the binary into the gravitational wave-dominated state before these torques align 
the binary with the angular momentum of the mass supply. 


Introduction 


The orbital evolution of binary black holes has been a subject o f considerable interest 
for more than thirty years, dating back to the work of Begelman et ah (|1980fh It has only 
grown in prominence in recent years, as it has become appare nt bot h that most reason ably- 


large galaxies contain supermassive black holes in their nuclei (JKormendy & Richstone 1995 
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Magorrian et al.l Il998l) and that today’s galaxies are the product of numerous mergers be- 
tween smaller galaxies ( Bell et aPl2006 : McIntosh et al. 2008 ; de Ravel et al. 2009 : Bridge et ah 


20101). each of which likely contains its own nuclear black hole. As pointed out bv iBegelman et al . 


(119801 ). when two massive black holes in a galaxy approach near enough one another to form 
a gravitationally-bound binary, interactions with external stars become much less effective 
at bringing them still closer together. A number of processes have since been suggested as 
viable mech anisms to push su ch a system closer toward the merger of the two black holes 
themselves ([Mer ritt fe M ilosavljcvic 20051 ). 


A particularly popular such mechanism ( 

Gould & Rix 

2000; 

Armitage & Nataraian 

2002: Escala et al. 2004. 

2005; 

MacFadven & Milosavlievic 2 

008l 

Cuadra et al. 

2009|) is ac- 


cretion of gas onto the binary black hole through a circumbinary disk. In principle, the 
orbital plane of such a disk could be wholly independent of the orbital plane of the binary, 
particularly at large radius. However, a number of pr ocesses have been poi nted ou t that 
tend to push the two orbital planes toward coincidence ([Nixon 2012; Miller & Krolik 2013). 

Whether the angular momenta of the two orbits align or counter-align, however, is less clear. 

Retrograde circumbin ary disks have certain special properties ( Nixon et al.lboil : Roedig fc Sesana 


2014: Ban kert et ah 2014). Because the angular momentum acquired by accretion is directed 
opposite to the binary orbital angular momentum, accretion in such a context tends to drive 
a rapid decrease in the orbital angular momentum while only modestly decreasing the or¬ 
bital energy. The result is an inc rease in the orbital eccentricity. Indeed, the first paper 
to consider evolution of this sort (INixon et ah 2 0111) concluded on the basis of qualitative 
arguments that this eccentricity growth causes the pericenter distance r p = a(l — e) to be¬ 
come so small that gravitational wave (GW) radiation rapidl y drives the system to merger. 
Although GW radiation generically decreases eccentricity (IPetersI 119641 ). in their view the 
binary orbit would retain non-zero eccentricity all the way to merger. 


Here we study this problem quantitatively. For a wide range of reasonable initial condi¬ 
tions and physical assumptions, we find the same robust result: rapid growth in eccentricity 
with constant apocenter a(l + e) over a time long enough to accrete a mass comparable 
to that of the secondary (for equal-mass systems accreting at the Eddington limit this is a 
Salpeter time, tg a i = 4.5 x 10' years, for a nominal radiative efficiency 77 = 0.1). At that 
point the eccentricity reaches e > 0.99, and gravitational wave evolution takes over. In this 
phase, GW losses occur almost entirely at pericenter, leading to evolution with constant 
a(l — e) until the orbit is nearly circular at a binary separation of ~ 100r 5 (the gravitational 
radius r g = GM/c 2 ). 

This evolutionary behavior can have profound effects on pulsar timing array (PTA) ob¬ 
servations of gravitational waves radiated by supermassive black hole mergers. If gas inflow 
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in the center of a post-merger galaxy takes place in a stochastic manner, retrograde disks 
should be roughly as common as prograde disks, but much more efficient at effecting the 
binary merger. Rapid evolution from initial separations of a ~ 10 5 r 9 down to a ~ 10 2 r n will 
great l y reduce the number of systems potentially observable in this regime (iKocsis fc Sesaiia 


2 0111; ISesana 2013 ; Rav i et a h 2014). As a result, the strength flDemorest et al.l l2013h of 
the stochastic GW signal in the PTA band / ~ 10~ 9 — 10” 8 Hz would be sharply di¬ 
minished. In addition, the early onset of GW-dominated evolution could lead to a decou¬ 
pling of the circumbinary disk from the inner binary at much greater radii than is typically 


i — —— ? _ — ~ — n _ ------" 

fMilosavlievic & Phinnev 

2005 

Chane et al 

2010; 

Tanaka & Menou 

2010; 

Noble et al. 

2012; 

Schnittman 

2013, 

2014 


"arris et al. 

2014bl 



Accretion-driven Evolution 


As shown in a companion paper c 


escribing ma gneto -hydrodynamic simulations of a ret¬ 


rograde circumbinary accretion disk ([B ankert et al. 2 0141) . the nearly axisymmetric density 


profile of the disk results in a very small gravitational torque on the binary. Thus, unlike 
the prograde case, the primary means of energy and angular momentum transfer from the 
disk to the binary is through direct accretion. 


The eccentricity e and semimajor axis a of the binary can be written as 


e = 


1 + 


and 


a = 


2EJ 2 (l + q) 6 
G 2 M 5 q 3 

q GM 2 


1/2 


(la) 


(lb) 


(1 + q) 2 2 E ’ 

with M the binary’s total mass, J angular momentum, E energy, and mass ratio q = 
M 2 /M l lei. The evolution equations are then given by 


e = 


2e 


Q ■ 


-2J/J - E/E + 5M/M + 3-- -q/q 


1 + q 


and 


1 — q 

a/a = —E/E + 2 M/M + -j— j —q/q. 


(2a) 


(2b) 


All of the terms on the right hand sides of both equations can be written in terms of the 
fractional accretion rate M/M. If the specific angular momentum of the accreting material 
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is jgas, then J / J = \j gaa /(J/M)]M/M. We find it convenient to parameterize j gaa in terms 
of the canonical Keplerian momentum via j gas = J\/GMa , giving 


j = 7 (i + g ) 2 i m 

J J q (1 — e 2 ) 1 / 2 M ' 


; From the MHD simulations of an equal-mass, circular retrograde binary in B a nke rt et ah 
(12014J) . we find J ~ —1.25. Assuming the inner edge of the circumbinary disk is limited 
by the radius of the binary orbit at apocenter, we believe it is reasonable to scale this term 
with eccentricity as J ~ —1.25\/1 + e. 


We parameterize the relative fraction of gas accreting onto each black hole as f\ and / 2 
(A + f 2 = 1), giving an expression for q: 


Q 

Q 


1 + q 


(/2 - qfi) 


M 

M 


(4) 


Here and in the following, quantities subscripted by 1 pertain to the primary, those sub¬ 
scripted by 2 pertain to the secondary. 

The relation between E and M is harder to determine because it depends on how much 
dissipation is associated with attaching the orbital streams to the members of the binary. 
Because shocks are extremely likely, dissipation will in most instances be very significant. 
Here we choose to estimate the energy delivered by accretion on the basis that the system is 
maximally dissipative, that is, that shocks are capable of dissipating as much orbital energy 
into heat (and likely radiating it) as is consistent with momentum conservation. I 11 other 
words, the gas streams inward from the circumbinary disk and then collides inelastically 
with one or the other of the “minidisks”, the individual accretion disks around each black 
hole. 


We assume that the accretion stream encounters the minidisk when the binary is at 
apocenter. Although we cannot confirm this assumption with our simulation data (the 
binary in the simulation had zero eccentricity), it is quite plausible: after all, this is the 
orbital phase at which the members of the binary are closest to the circumbinary disk, and 
also the phase at which the binary spends most of its time. 

At apocenter, the total energy in the orbit can be written as 


E 


GM X M 2 GM 2 q 


2a (1 + qf 


2 a 

GM 2 q 1 2 1 2 

a(l + e) (1 + q) 2 2 1 1 2 2 ’ 


( 5 ) 











5 


where in the second line we separate the energy into potential and kinetic terms, with i>i i2 
the velocities at apocenter, and r 1)2 = a(l + e)(q, 1)/(1 + q ) the radial positions of the two 
black holes. Here we adopt a notation where the respective quantities for the two black holes 
are written (1, 2). The rate of change in the binary’s potential energy is then 


E p = 


G{M 2 M l + M 2 Mx 
a(l + e) 


GM 2 M 


q 


fi + 


a(l + e)M\l + q 1 + q 


-/ 2 


= 2 E 


M 1 1 + q 

M 1 + e q 


( 9/1 + h) 


( 6 ) 


For the kinetic energy terms, it is convenient to write the kinetic energy at apocenter 
in terms of the linear momentum of each BH: E\~ 1,2 = l/2(p 2 2 /Mi ) 2 ). Then we have 


Tp _ Pl,2Pl,2 Pi,2 , . 1 2 71 '/ff 

Eh 1 5 — “TV- _ , Mfi t 2 — Vi t 2Pl t 2 — 2 V 1,2™J1,2 


vfc1 ’ 2 Af 1)2 2M 2 2 

where the linear velocities at apocenter are 


GM 1 — e (q, 1) (q, 1) 

Vi,2 = \ —r\T— T—„ = ^po- 


a Vl+el+g 


1 + q 


(7) 


( 8 ) 


The momentum change comes from the accretion stream: 

Pi,2 — %isl,2 Mfip 

which has velocity 


%ssl,2 — 


Jgas 

r l,2 


(9) 


( 10 ) 


As above, j gas is the specific angular momentum of the gas in the accretion stream. We can 
now write 


3 V GMa 1 + q I GM jl — e I 1 1 + q I 1 1 + q , , 

W-1.2 = , x , ,1 = j\ —— \ 7TTT = 3 v,po\l- ^7—IT • ( n ) 


a(l + e) (q, 1) 


a V 1 + e V 1 — e 2 (q, 1) 


1 - e 2 (q, 1) 


Now we can find the change in kinetic energy: 


E k l,2 = v. 


apo \/ 1 _ e 2 
= —Mfi 2 v. 




/ 

apo 




l-e 2 2(1 + 9)= 


( 12 ) 
































Writing 


(13) 


V — 
^apo 


—2 E\ 1 — e (1 + qY 


M J 1 + e q 

and combining this with the expression for the potential energy, we get 


-Ei ,2 M „ 1 

IT = M h ' 2 1 + e 


-2 J 


l-e(l + g) 2 1 + q, 

T—-+ (1-e ){q,q + 2- (q,l) 

1 + e q q 


(14) 


Recall that, since J < 0, this expression is always strictly positive. Therefore the binary 
always evolves toward greater binding energy. 

We can now inspect equations (I2a|2bjl in a few limiting cases: g = 1, g C 1, e € 1, 
1 — e = 5 -C 1. The leading-order terms in each case are shown in Table [I] For q - Cl 
we assume f 2 = 1. Not surprisingly, the eccentricity always grows and the semi-major axis 
always shrinks for binaries in the presence of retrograde accretion disks. Also not surprisingly, 
the effect is fastest when q is small. Note that all evolution terms are proportional to 
q ^ 1 M/M, so the characteristic time scales can be rather shorter than the accretion time 
M/M = 4.5 x lCErh -1 yr, with m the Eddington-scaled accretion rate for a radiative efficiency 
of 7 j — 0.1. 


Table 1: Leading-order terms for e and a for limiting values of q and e. 


Q 

e 

e 

a/a 

1 

< 1 

1S— 

16 M 

—1S— 

16 M 

< 1 

< 1 

9 M 

2 q M 

9 M 

2 q M 

1 

1-6 

i oVs% 


< 1 

1-6 




Because we assume that all accretion takes place at apocenter, and the accretion flow 
is tangent to the binary orbit at that point, the impulse approximation implies that a(l + e) 
should remain constant. When e -C 1, this implies that e = —a/a, while e = — 2d/a when 
e ~ 1. This results are consistent with the limiting behavior shown in Table [I] 


3. Gravitational Wave-Driven Evolution 


In the previou s section, we saw that the binary eccentricity can grow rapidly via ret¬ 
rograde accretion. Nixon et al. 120111) also reached this conclusion, and predicted that, as 
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eccentricity grew and the pericenter separation decreased, gravitational waves would take 
over as the driving evolution ary me chanism of the system. Here we present a more de¬ 
tailed analysis of this process. Peters! ( 1964 ) gives the gravitational-wave-driven evolution of 


eccentricity and separation in the post-Newtonian regime: 


(a/a) 


64 


GW 


Q 


(ejcw = 


5 (1 + q) 2 
304 q 
' 15 (1 + q) 


(1 


= 2n-7/2 


73 


37 


1 + 777^ + — e 


24 


96 


.-i 


.(1 _ e 2 )- 5 / 2 e 


1 + 


121 

304 6 


^ t 


-l 


(15a) 

(15b) 


We have scaled the spatial dimensions by the gravitational radius r g = GM/c 2 and the time 
scale by t g = r g /c. 



a(r fl ) 


Fig. 1.— Characteristic time scales for various evolutionary processes. As the system 
evolves from right to left, first the eccentricity rises rapidly with constant a(l + e) and 
then gravitational-wave-driven evolution takes over, decreasing both a and e, while keeping 
a(l — e) roughly constant. We take for initial values M = 10 8 M 0 , q — 1, m = 1, a 0 = 10 5 r s , 
and e 0 = 0. t e<gw = f a;gw over a broad range of semi-major axes because a(l — e) is constant. 


For a given a, e we can now define the characteristic time scales t e = (1 — e)/e and 
t a — a/a for the competing evolutionary processes: accretion and gravitational waves. For 





















a fiducial binary system with initial parameters M = 1O 8 M 0 , q — 1, fn — 1, a = 10 5 r g , and 
e = 0, we plot in Figure [[] the relevant timescales as a function of a(t). The graph should be 
read from right to left, as a shrinks with time. Whichever timescale is shortest at a given 
orbital separation is the dominant process governing the binary evolution at that point. 

Consider, for example, the evolution of our fiducial case. As it accretes from its cir- 
cumbinary disk, its eccentricity grows rapidly for a few million years, while its semi-major 
axis simultaneously shrinks by a factor of two. At this point, the high eccentricity leads 
to an extremely sharp decrease in t Qjgw and the system transitions from accretion-driven to 
gravitational wave-driven evolution. For highly eccentric orbits, most of the gravitational 
radiation is emitted at pericenter passage, when the two bodies have the highest velocity. 
This fact justifies an impulse approximation in which the energy and angular momentum 
are lost instantaneously at pericenter, so the binary evolves with constant a(l — e). 



1-e 


Fig. 2.— Evolution of a and e for the fiducial binary model with M = 1O 8 M 0 , q — 1, and 
fn — 1 (solid curve). We also show variations from this model with M = 10 9 M 0 (dotted 
curve), q — 0.1 (dashed curve), and fn = 0.1 (dot-dashed curve). 


In Figure [2] we show the mutual evolution of a and e for a few variations of the fiducial 
model parameters. In all cases the binary begins in the upper-right corner with a = 10 5 r g and 
e = 0. The dependence on the various parameters is straight-forward. Because t g grows with 
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M, whereas for fixed rh the accretion time is mass-independent, more massive systems reach 
a higher eccentricity before gravitational waves begin to dominate. For decreasing rh, the 
opposite happens: the accretion time increases even though the gravitational time remains 
constant, so the system leaves the accretion-dominated regime earlier, at a smaller eccen¬ 
tricity. For small q, the gravitational wave timescales increase and the accretion timescales 
decrease; both effects combine to push the system to higher eccentricity before entering the 
gravitational wave regime. More generally, because a(l + e) is constant so long as the evolu¬ 
tion is accretion-dominated, and the transition to GW-driven evolution occurs when e ~ 1, 
the semimajor axis at this transition is always ~ 1/2 the initial semimajor axis. 

We can solve for the peak eccentricity e* by simply solving e acc = —ecw hi the limit of 
1 — e* = i < 1. For e acc we use equations (I2a|) and (ED and the observation that the J/J 
terms dominate the evolution at high eccentricity. In this limit we find 


<5* = 4 x 10“ 3 


.(1 + 9) 


m 4 M a 1 rh 1 


4^5 


1/3 


(16) 


which agrees quite well with the numerical evolutions plotted in Figure [2j Here we use the 
scaling expressions a 5 = a/(10 5 r g ) and M 8 = M/(10 8 M Q ). For a wide range of parameters 
we find the time f t rans needed to evolve from e 0 to e* can be approximated as 


t 


trans 


1 _ Q_ 

3 (1 + q) 2 


(1-eo ) 1/2 m % a \. 


(17) 


Because q/( 1 + q ) 2 , also called “the symmetric mass ratio”, is always < 1/4, ftrans < 
(l/12)(fs a i/rh), where (tsai/rh) is the time over which accretion increases the mass of the 
system by one e-fold. 


4. Gravitational Wave Signal 


For binaries on circular orbits, the gravitational waveform is a sinusoidal signal with 
frequency exactly twice that of the orbital frequency / or b- For non-zero eccentricities, the 
waveform is decidedly non-sinusoidal, as the instantaneous orbital frequency increases as the 
black holes move from apocenter to pericenter, peaking at 


fp ./orb 


(1 + e) 1 / 2 
(1 - e) 3 / 2 


(18) 


Consequentl y, the range of harmonics carrying substantial power in creases with increasing 
eccentricity (ITessmer fe Gopakumarl 120071: lEnoki fe Naga,shirnall2007lh 
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We calculate the amplitude h(t) of the leading-order GW waveform from the time- 
varying quadrupole moment: h(t) ~ Q. The power emitted in the waveform is Tcw(f) ~ 
h 2 (t). Taking the Fourier transform Lgw(/)i we get the distribution of power in the indi¬ 
vidual harmonics. To get the correct normalization for Lqw(/), we use the orbit-averaged 
energy radiated dE/dt 

dE _ GM 2 q d 

IF ~ (1 + q) 2 ~2a? ’ ^ ' 


where a/a is given by [Peters] (1l964h . We hnd that, for 1 — e -C 1, the instantaneous dE/df 
peaks around f p . For circular orbits /gw = 2 f p , so the fact that /gw ~ fp for eccentric 
orbits reflects the fact that the GW radiation is emitted over a finite part of the orbit 
around pericenter, during which / < f p . 



f (Hz) 


Fig. 3.— Gravitational-wave power spectrum for the fiducial binary model with M = 1O 8 M 0 , 
q — 1, and m = 1 (solid curve). We also show variations from this model with M = 1O 9 M 0 
(dotted curve), q — 0.1 (dashed curve), and m — 0.1 (dot-dashed curve). The thin dashed 
line is dE/df ~ / -1//3 . 


Figure [3] shows the GW power spectrum for the same four cases considered in Figure 
[2} fiducial model with M = 10 8 M 0 , q — 1, and m — 1 (solid curve), along with variations 
from this model with M = 1O 9 M 0 (dotted curve), q = 0.1 (dashed curve), and m = 0.1 
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(dot-dashed curve). The thin dashed line corresponds to dE/df ~ / 1//3 , the classical result 
for purely circular evolution. 

Figure [3] can best be understood by close comparison with Figure [21 As the binary’s 
eccentricity increases through accretion, the GW signal rises sharply in amplitude and fre¬ 
quency until the binary evolution switches over to being dominated by gravitational waves. 
At this point, the binary evolves with constant a(l — e), and thus constant /gw, building up 
power in a single peak. Eventually, the eccentricity is small enough that the system follows 
the standard circular evolution track, with both amplitude and frequency evolving to give 

dE/df ~ r 1 / 3 . 

Comparing it with Figured! we can also easily understand the relative offsets between 
the different models plotted in Figure [3J For the more massive system with M = 1O 9 M 0 
(dotted curve), the peak frequency shifts to the left since the orbital frequency decreases with 
increasing mass. At the same time, the GW amplitude increases because dE/df oc M 5//3 . For 
smaller accretion rate (m = 0.1; dot-dashed curve), the peak eccentricity is smaller, giving a 
lower peak frequency [see eqn. flT%D ], but it eventually evolves along exactly the same curve 
as the fiducial case. For smaller mass ratio (q = 0.1; dashed curve), the peak eccentricity is 
larger, leading to a higher peak frequency. The amplitude scales like dE/df oc q, so the net 
effect is to shift the fiducial curve down and to the right. 


The area under the solid black curve is actually the same as that under the thin dashed 
line because they both represent the total energy radiated going from a ~ 10 5 r 9 down to 
merger. Thus, all the power normally radiated at low frequencies is concentrated in a small 
peak around / ~ 10“' — 10“ 6 Hz. This modification of the GW spectrum can potentially 
have a major impact on the stochastic gravitational wave signal being searched for by pulsar 
timing arrays (PTAs), which are most sensitive to GW signals in the 10“ 9 — 10“ 8 Hz range 
( Hobbs et al. 2010 : Shannon et ah 2013 : Demorest et ah 2013 : Sesana 2013). 


If we take the initial separation to be somewhat smaller (a = 10 4 r g ), we find qualitatively 
similar behavior in the evolution, but from equation (TT6|) we see the peak eccentricity is 
“only” e ~ 0.9 — 0.95. In this case, the peaks in Figure [3] all shift to lower frequencies by 
a factor of ~ 2. Integrating over the cosmic population of SMBHBs, this shift could give 
a PTA stochastic signal with a peak around / ~ 10“' Hz and an amplitude roughly 50% 
higher than the canonical value, while sharply reducing the signal at lower frequencies. 
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5. Rapid Approach to Decoupling 


During the accretion-driven phase, the binary separation a(l + e) remains constant, 
so the radius at which the secondary is farthest from the center-of- mass, and therefor e best 
able to capture accre t ing gas, remains unchanged. Numerous studies (IRoedig fe Sesanall2014l: 
D’Orazio et al.ll2013t iFarris et al.ll2Q14bl: IShi fe Kroliki 12 01 5h have found that so long as q is 
not too small, the secondary captures the majority of the accretion flow from a surrounding 
disk in prograde binaries; that fact is likely to carry over to the retrograde case. On the other 
hand, during the GW-driven phase, r p = a(l — e) is fixed, while the apocenter evolves like 
v a = r p ( 1 + e)/(l — e), which falls from ~ 2a* to ~ a*(l — e*) as the orbit grows less eccentric 
(a* and e* are the semi-major axis and eccentricity at the point when the GW-driven phase 
begins). 

The rapid diminution in apocenter distance during the GW-driven phase may be faster 
than the characteristic inflow rate through the inner region of the accretion disk. If during 
the relevant phases of orbital evolution, the disk is in the Shakura-Sunyaev “outer zone” (gas 
pressure-dominated, free-free opacity greater than electron scattering opacity), this inflow 
time at radius r = rs/(10 5 r 9 ) is 

7 x 10 8 m- 3 / lo alt /5 M 8 6/5 r 5 7/5 yr. 


^inflow 


( 20 ) 


Here a = 10 _1 ck_i is the parameter defined by IShakura fc Sunvaevl (119731 ) to describe the 
ratio of the vertically-integrated internal stress to the vertically-integrated pressure, scaled 
to its canonical value of 0.1. While f in fl ow is somewhat larger than both t aACC and t e ^ cc during 
the accretion-dominated phase (see Fig.[T]), this does not actually imply decoupling, since the 
important quantity for comparison is the apocenter evolution timescale, which is essentially 
infinite at this phase in the evolution, so the disk is likely in inflow equilibrium. 


Equation (1201) assumes a steady-state thin disk geometry, in which the local temperature 
is controlled by a balance between heating due to accretion and local radiation. However, it 
has been recognized for ma ny years that for typical AGN parameters, such disks are unstable 
to self-gravity at r > 10 3 r „ flShlosman fe Begehnanlll987l : lGoodrnanll2003l : I Jiang fc Goodman 


2011; Dunhill et al. 


2 0141) . and this range of radii is of particular interest to the questions 


studied here. Unfortunately, little is known about what actually happens as a result. It is 
possible, as these papers have discussed, that the disk fragments and a significant part of its 
mass is concentrated into new stars. 

On the other hand, we have observational evidence at larger scales (> 10 5 r p ) that ac¬ 
creting material forms geometrically thick structures, so-called “obscuring tori;” most recent 
work has supposed th at they are supported by the AGN's radiation pressure interacting with 
dust grains in the gas ( Pier fe Krolik 1992 ; Chang et alil2007 : Krolik 2007 ; Dorodnitsvn et al. 
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20111. 2012; Roth et ah 2012 '). To account for additional mechanisms that might alter the 
disk’s vertical support, it is useful to rephrase the inflow time in terms of the disk aspect 
ratio h/r, its ratio of vertical scale height to radial position. The inflow time then becomes 


inflow = 5 x 10' a_\M s r/ 


„3/2 


h/r 

OOl 


-2 


yr. 


( 21 ) 


We scale h/r in this instance to 10 2 because the aspect ratio predicted by classical disk 
theory is 2 x 10” 3 CK _1//10 m 3//20 M 8 _1 ^ 4 r^ 20 , and this is presumably a lower bound. 

By comparison, the characteristic gravitational wave timescale tew — a/a in the limit 
of e —» 1 is 

(i + q) 2 - - * ( x 7/2 


tew — 5 x 10 5 -- ——Mg, at 

q 5 V 0.003 


yr. 


( 22 ) 


We have scaled to the typical orbital elements when gravitational wave evolution begins to 
dominate. We therefore expect that in most instances, the apocenter changes too rapidly 
for the bulk of the disk to follow once GW emission controls the system’s evolution. 


Milosavlievic & Phinncy (120051) argued that whenever the binary evolution does become 


faster than the disk’s inflow time, mass inflow is cut off, “decoupling” the disk from the binary 
and ending electroma gnet ic emission due to accretion. On the other hand, the global MHD 
simulation of [Noble e t al. (1201 2) showed that the accretion rate can, in fact, be sustained 
even in this regime—edge effects can substantially raise the ratio of stress to pressure above 
its usual value. Similar effects are gen erically seen during the transient phases of almost 
all global MHD accretion simulations (lHawlev et al. 201 ll; ISorath i a et al.l 2012). In fact, 
analogous edge effects occur even in purely hydrodynamic treatments flFarris et al.ll2014al ). 
Thus, the degree to which mass accretion actually diminishes when the orbital evolution 
time is shorter than the disk inflow time is uncertain. In either case, retrograde disks should 
reach this “decoupling” considerably sooner than prograde disks. 


6. Disk-Binary Alignment 


A binary and its disk can remain in a stable counter-aligned orientation if Jd is k < 2 Jb; n 
(Ni xon 12012). However, this criterion is difficult to satisfy in the long-run. The relevant 
disk mass is the total amount of mass passing through it integrated over time, and, as the 
binary’s eccentricity approaches unity, Jbin falls. Both effects make it easier to f lip the binary 
around while still conserving the total angular momentum of the system (iRoedig fe Sesana 
2 0141 ). We have shown already that the angular momentum conveyed with mass transfer 
is insufficient to accomplish this transition to a prograde orientation before gravitational 
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wave emission brings the binary to merger; however, the situation changes when the binary 
orbital plane and the disk plane are not exactly co-planar because in that configuration 
angular momentum can also be transferred between the two systems via the mutual torque 
arising from the binary’s quadrupole moment. 


Suppose, then, that the angular momentum fed to the disk at large radius is aligned 
not exactly opposite to t he ang ular momentum of the binary, but only approximately. As 
shown by Miller & Krolik (1201311 . when the binary orbit is circular, the inner portion of the 
disk aligns relatively quickly with the orbital plane, but beyond a certain distance the disk 
maintains its intrinsic orientation. Once the disk has reached a steady-state shape, the rate 
at which angular momentum is delivered from the disk to the binary is ~ Mjr, where jr is 
the magnitude of the misaligned specific angular momentum in the disk at the radius where 
the disk orientation transitions from aligned to misaligned, and the radius of that transition 
is Rt ~ (r/h)ra, where ( h/r)x is the disk aspect ratio at Rt- 


For the present application, this result must be generalized in two ways. The first is that 
their work assu med the binary has a circular orbit, so it must be extended to eccentric orbits. 
Following Landau fe Li fshitzj (I19691) . it is easy to show that the time-averaged quadrupole 
moments of an eccentric binary are 


Qxx 

Qyy = 


M 


q 


2 (1 + g ) 2 

M q 

(l + <?) 2 


a 2 (l + 4e 2 ) 
a 2 (l — e 2 ), 


(23a) 

(23b) 


where we define the x— and y— directions as parallel to the major and minor axes of 
the ellipse, respectively. If the disk is tilted in the e — z plane, where z _L x, y and e = 
cos <f>x + sin (fry , the torque is proportional to Q xx cos 2 </> + Q yy sin 2 </>. 


Because the quadrupolar torque is precessional, both the orientation of the disk tilt 
plane </> and the misalignment angle 9 can be functions of radius within the disk. The total 
binary-disk torque is then 


T = 


37 r q 

"2" (i + q ) 2 


fi 2 a 2 sin 9 0 R^T t x 


dxx 


_ 2 sin 9(x) T(x) 
sin On Sr 


[(1 + 4e 2 ) cos 2 4>(x) + (1 — e 2 ) sin 2 4>(x)\ , (24) 


where quantities with subscript T are evaluated at Rt, beyond which the disk’s intrinsic 
alignment makes an angle 6q with the binary. The orbital frequency at this location is Qt, 
Tit the surface density, and x = r/R t . The form of this expression for T is wr itten so 
as to make the integral both dimensionless and order unity, as in Miller fe Krolik (20131). 
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Following their notation, we designate it as Z, but now recognize that, unlike in the case of 
circular binary orbits, it is a function of e and precession phase profile </>(r). 


The second generalization is a consequence of the implicit assumption in Miller fe Krolik 


f l2013n that alignment of the binary remains incomplete even after the alignment front has 
propagated all the way through the inner disk, so that the torque is concentrated near 
r ~ Rt . If, when the system is first formed, the misaligned angular momentum contained 
in the disk within a radius smaller than Rt (the steady-state alignment transition radius) is 
comparable to the angular momentum of the binary, the binary will have reoriented before 
the alignment front arrives at the steady-state radius. 


In fact, just this situation may often occur in the context of binary black holes. Sup¬ 
pose that the circumbinary disk extends from an inner radius ~ a (because it is obliquely 
retrograde) out to an outer radius r out . As we did in the previous section when discussing 
the approach to decoupling, we will describe the disk’s surface density profile in terms of a 
particular aspect ratio h/r, rather than calculate the aspect ratio from a model for the disk’s 
thermal state. The ratio between its misaligned angular momentum and the binary’s total 
angular momentum is then 


Xdisk „ ■ a ™ ( / ',3/2/ /l \2 [(^out/ a ) 1](1 d ?) G r g 

—— = 27Tsm0 o — a r g y 1 r h) - - - —m - 2 

Jb in q( 1 - e 2 ) 1 / 2 n T c 2 


(25) 


where kt is the Thomson opacity. If the scale height varies as a function of radius, ( r/h) 2 
should be regarded as the angular momentum-weighted mean of (r/h) 2 . Evaluated for our 
fiducial values, this ratio is 


Jd i 


disk 


Jbin 


~ 0.06 sin 0 o 






1 (1 + q ) 2 

(1 — e 2 ) 1 / 2 q 


[(r 0 nt/a ) 2 - l] . 


(26) 

In other words, for our fiducial values, the disk does not need to extend to a large multiple 
of a in order for its angular momentum content to match the binary’s. However, it is also 
clear that this ratio is sensitive to a variety of parameters. 


Suppose first that Jdisk ^ Jbm- In this case, as the alignment front propagates out 
through the disk, the corresponding torque on the binary causes its angular momentum to 
swing toward the direction of the disk’s angular momentum. Because the direction of the 
binary’s angular momentum is initially nearly opposite to that of the outer disk, this swing 
initially increases the angle between the two orb ital planes, the inclination angle growing 
exponentially if it is initially small f Scheuer fc Feiler 1996 1^1. 


1 Nixon (2012|) claimed otherwise, but he assumed a disk without continuing accretion supply and ran for 

too short a time to test stability. 


















- 16 - 


However, equation (126|) suggests that there could also be disks of limited radial extent 
containing less misaligned angular momentum than the binary, and the maximum r out /a for 
which J disk < J bin depends on the accretion rate, binary mass, etc. When r out is small enough 
for J disk < J bin , the disk counter-aligns swiftly, i.e., it moves into the binar y orb ital pl ane. 


b ut o rbiting retrograde to the binary motion. Adapting the formalism of iSorathia et al . 


1 20131) from the Lense-Thirring problem they considered to the closely analogous problem 
of alignment with a binary, we find that the alignment front moves outward at a rate 


v f {r) ~ 


3 q a 2 fl(r) 
8 (1 + q) 2 r 


Z(r) 


(27) 


In this expression X is redefined so that it covers the region within r rather than within Rt- 
Consequently x in equation ()24l) becomes r'/r and Tit becomes T(r). Thus, for order unity 
mass ratio in the binary, Vf ~ (a/r) 2 u orb (r), not much slower than the orbital speed, and the 
time for alignment to be achieved throughout a truncated disk is ~ (r out /a) 2 binary orbital 
periods at its outer radius. 


Once such a disk is entirely aligned, the binary exerts no further torque on the disk 
proper, but does continue to interact with material at and beyond its outer boundary, r out . 
Because the torque decreases rapidly with increasing radius (the precession frequency is oc 
r - '/ 2 ), only the nearest matter is significant. Typically the mass just outside the disk should 
be much smaller than the mass within the disk because the surface density of the disk is ~ 
200 (a/ 0 . 1 ) -1 (m/r 7 )[(r/h)/ 100 ] 2 (r/ 10 a)“ 1 / 2 (a/ 10 5 r 9 ) -1 / 2 /v^ 1 , a surface density much greater 
than would be expected from the surrounding interstellar medium. Consequently, the torque 
that the binary continues to exert is on the material that has just become attached to the 
disk, arriving with the intrinsic orientation of the mass supply: T = sin doMr 2 ut f](r out ). 
This torque is smaller than the steady-state torque estimated through the formalism of 


Miller & Krolik (2013j) by the ratio (r out / Rt) 1 ^ 2 . 


Whether the binary evolves to gravitational wave-driven merger faster than it aligns 
with the external mass supply depends on whether t a ii g n, the time required for this torque to 
change the orientation of the binary, is larger or smaller than Arans, the time for accretion to 
drive binary orbital evolution. Because both timescales depend on the angular momentum 
content of the binary and the mass accretion rate, their ratio is very simple: 

= 1 sin 6^ (r out /a) 1/2 . (28) 

Olign O 

Thus, the maximum size of the circumbinary disk permitting retrograde orbital evolution to 
reach its conclusion is limited by two constraints. One is that the timescale for eccentricity 
growth must be shorter than that for disk-binary alignment: 


rout < (9/ sill 2 0 o )a. 


(29) 
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For moderate obliquity, r OVLt /a could be as large as ~ 10-100 and still satisfy this condition. 
The other constraint is the one described above in equation ([26]) . namely that Jdisk be < Jbin- 
This requirement depends on more parameters. For our fiducial choices, it gives a tighter 
constraint, r out < 2a/-\/sin 9 0 , but this constraint becomes similar to or looser than the 
timescale condition when ( m/r])M 8 (a/10 5 r g ) 3 ^ 2 [(r/h)/100] 2 (l + q) 2 /q < 0.1 (e.g., smaller 
accretion rates, more compact binaries, or thicker outer disks). 


ft is worth pointing out that the timescale criterion we h ave us e d can be rephrased in 
terms of the total angular momentum and mass budgets. K ing et ah (120051) argued that the 
binary can be flipped when Jdisk/Aiin > 2. This criterion applies when there is no resupply 
of mass to the disk. The alignment timescale in our case should therefore be posed in terms 
of the ratio of the binary angular momentum to the rate at which (misaligned) angular 
momentum is supplied to the disk. The transition time t trans is defined above in equation 
cm, and is roughly the time it takes to accrete 20% of the secondary’s mass. 


Up to this point, we have ignored the question of whether such truncated disks might 
be expected in Nature. In fact, they may be reasonably common. Presumably what defines 
the outer edge of the circumbinary disk is the point at which interstellar gas, following 
ballistic orbits in the gravitational potential created by the binary black hole and the host 
galaxy, shocks against the disk’s outer edge and joins the accretion flow. Another reason 
why the inventory of mass immediately outside the disk should be much smaller than the 
disk mass is, in fact, that mass in this region travels inward ballistically, whereas within the 
disk inward motion is dependent upon the comparatively slow process of internal angular 
momentum transport. The place where external matter joins the disk could be at a radius 
as small as that for which the specific angular momentum of the interstellar gas supports a 
roughly circular orbit; alternatively, it could be at larger radius if the accretion flow (having 
this orientation) has lasted long enough for the outer edge of the disk to be pushed outward 
significantly by the outward angular momentum flux of accretion. 


Details of the formation of binary black holes are complex and depend upon numerous 
parameters of the galaxy mergers preceding the formation of such a binary. However, chaotic 
motion of interstellar gas on the scale of parsecs to tens of parsecs with orbital velocities 


o f sev eral hundred km/s can often be seen in simulations of specific scenarios ( Maver et al. 


20071 : [Roskar et al. 2014). By comparison, the physical scale of the of the disks we have 
been considering is similar or somewhat smaller, ~ 5 pc for our fiducial parameters with 
Amt ~ 10a. Thus, a truncated obliquely retrograde disk of the sort we have discussed may 
well be formed. 
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7. Discussion 


To form a circumbinary accretion disk, one requires a SMBH binary and a supply of gas. 
Galactic mergers are almost certainly required to provide the binary, but although mergers 
may enhance the availability of gas, they are not strictly necessary, as we see plenty of 
examples of AGN in non-merging galaxies. Assuming a merger did take place and a SMBH 
binary has formed, if subsequent nuclear activity is driven by the stochastic accretion of 
molecular clouds with random orientation relative to the central SMBHB orbit, we expect it 
would be equally likely for circumbinary accretion disks to form in the prograde or retrograde 
orientations. 


Wi th major mergers happening every ~ 3 — 10 Gyr per galaxy in the local universe 
flLotz et al.l 20 111), a nd ty pical AGN lifetimes ~ 10' years flMart ini ,20041) with duty cycle ~ 
ICG 2 flMarconi et al.l 20041). we should expect multiple accretion phases between mergers. For 
prograde circumbinary disk s, doub ling: the secondary mass leads to the binary orbit shrinking 
by a factor of order unity (jShi et al. [2012)), so many prograde episodes would be necessary 
to reach the GW-dominated phase. However, retrograde disks lead to much more rapid 
orbital evolution: a mere 20% increase in the mass of the secondary suffices. For the same 
accretion rate at large radius, retrograde disks are more efficient in driving binary orbital 
evolution than prograde disks not because they tran smit more mass to the binary [both cas es 
maintain accretion flows close to inflow equilibrium (iFarris et al.ll2014bt IShi fe Krolikll2015f) ]. 
but because retrograde accretion changes the specific angular momentum of the binary by 
a much larger factor. The angular momentum carried to the binary by accretion in the 
prograde case can in fact significantly counteract the an gular m omentum lost by the binary 
through gravitational torques on the circumbinary disk flSh i et al.l 120121. 


Thus, within this paradigm of multiple, stochastic accretion episodes between mergers, 
there is a strong statistical bias toward such rapid eccentricity evolution via accretion from 
retrograde disks that the GW-driven state is achieved quickly. When the mass accreted 
during a single episode of constant angular momentum orientation is less than roughly half 
the secondary mass, substantial orbital evolution occurs only for roughly retrograde orien¬ 
tation. Moreover, when substantial orbital evolution does take place due to a retrograde 
accretion episode, it drives the binary toward merger, so that no further episodes, retrograde 
or prograde, occur. 


Torques couple binaries to obliquely-oriented circumbinary disks. When J disk ■ ^binary < 0 
and enough mass is available (either in the initial disk or over the course of its feeding) that 
the total disk angular momentum is larger than the binary’s, these torques can “fl ip” the bi¬ 
nary orbit so that it_aligns with the disk plane and orbits in the same sense OScheuer fe Feiler 


19961 : King et al. 120051 1. However, we have shown that if the disk is fed at radii not too much 
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larger than the binary semimajor axis, the angular momentum instantaneously stored in 
the disk can be small enough, and the evolution of the binary orbit caused by retrograde 
accretion can be rapid enough, that the binary is driven to the gravitational wave-dominated 
stage before its orbit is flipped. Thus, the condition for retrograde accretion to cause merger 
before the binary orbit is reoriented amounts to a condition on the disk’s outer radius relative 
to the scale of the binary orbit. 

As a result, high-eccentricity and ultimately GW-driven SMBHBs may be a common 
result of gas-driven black hole mergers unless any of three conditions applies: accretion 
episodes frequently occur in which the total mass accreted is large enough to cause substantial 
binary orbital evolution even if the rotation sense is prograde; the orientation of the accreted 
matter’s angular momentum is consistently correlated with that of the binary black hole 
system; or the outer edge of the circumbinary disk is so far outside the binary that retrograde 
systems are flipped before their orbits evolve. In evaluating this statement, we emphasize 
the qualifier “gas-driven” because gas accretion may not always be the principal element in 
binary black hole evolution: stellar interactions may be more important in some cases. This 
is especially relevant for the massive gas-poor elliptical galaxies that host the most massive 
black holes, and which are thought to dominate the PTA signal. These galaxies also exhibit 
decreased nuclear activity. It is also possible that even when stellar interactions are weak, 
the evolution of the circumbinary disk at large radii may be complicated by self-gravity and 
internal star formation. 


One consequence of the gas-driven, rapid merger scenario would be a distinct absence 
of binary quasars with separations between ~ 0.01 — 1 p c, a source population that has 
been notoriously difficult to detect flEracleous et all 2 01 2). To the degree that gas-driven 
accretion dominates binary evolution in PTA sources, the effects we have pointed out would 
have profound effects on the potentially-observa ble GW signa l, effectively moving most of 
the GW power from the nHz to the /iHz band flSesana 2013). By the time such systems 
reach the mHz band, they should be highly circularized. 


Because the GW evolution time drops so rapidly with increasing eccentricity, it is likely 
that at some point before black hole merger, the circumbinary disk will no longer be able to 
stretch inward as rapidly as the binary apocenter distance shrinks; that is, the the system 
will “decouple.” Yet even in such a situation, there still might be enough low-density, 
highly magnetized gas in the central gap to lead to an observable EM signal. Understanding 
the precise nature of such a system will require new 3D global MHD simulations that can 
sufficiently resolve both the disk and the gas flow onto the individual black holes over long 
periods of time. 
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